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Two distinct phase transitions occur at different temperatures in QCD with adjoint fermions 
(aQCD): deconfinement and chiral symmetry restoration. In this model, quarks do no explicitely 
break the center Z(3) symmetry and therefore the Polyakov loop is a good order parameter for the 
deconfinement transition. We study monopole condensation by inspecting the expectation value 
of an operator which creates a monopole. Such a quantity is expected to be an order parameter 
for the deconfinement transition as in the case of fundamental fermions. 
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1. Introduction 

As lattice simulations seem to show, chiral symmetry restoration and deconfinement occur for 
QCD at the same temperature, making it extremely hard to disentangle them. On the other hand, 
QCD with quarks in the adjoint representation of SU (3) (aQCD) is a model in which two transi- 
tions take place at different temperatures In aQCD, the Z(3) center symmetry is not broken by 
adjoint fermions and therefore the Polyakov loop is an order parameter even at finite quark mass as 
in the pure gauge case. 

The authors of ref.[jl|] find on the lattice two distinct phase transitions; they observe a first order 
deconfinement transition and a continuous chiral transition. They also check that the observables 
which are sensitive to deconfinement are not significantly affected by the chiral transition. 

Another way to study confinement consists in looking for magnetic charge condensation, 
which signals dual superconductivity of the vacuum [Q, ||, |], ^]. In this case, one constructs an 
operator which carries magnetic charge and looks at its vacuum expectation value, which is ex- 
pected to be different from zero in the confined phase and strictly zero in the deconfined phase. 
Our goal is to study monopole condensation in aQCD in order to check that it is unaffected by the 
chiral transition. 

In this report we present the results of a lattice simulation of aQCD in which we started our in- 
vestigation; it has to be considered as work in progress. The qualitative behaviour of the monopole 
order parameter was found to be consistent with the expectations; however, for time limitations, 
we could not simulate quark masses that were light enough to study its chiral properties. We are 
currently continuing our investigation in this direction. We also need to understand the scaling 
properties of the monopole condensation parameter. 

1.1 aQCD 

Quarks in the adjoint representation of QCD have 8 colour degrees of freedom and can be 
described with 3x3 hermitian traceless matrices 

Q(x) = Q a {x)X a (1.1) 

using Gell-Mann's X matrices. In order to write the fermionic part of action for this model, the 
8— dimensional representation of the gauge links (which is real) must be used: 

Ufa = ±Tr(X a U i3) X b ul 3) ) (1.2) 

The full action is therefore given by 

S = S G [U {3) ] + 1 £Q(x)M(U {8) ) XJ Q(y) (1.3) 

x,y 

where So is the usual SU (3) gauge action and M is the fermionic matrix. The Polyakov loop is 
defined as in the SU (3) Yang Mills case 

L 3^(^HTrn[/ ( S 3 W)|> d- 4 ) 

x *o=l 
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This quantity is an order parameter for the spontaneous breaking of the center symmetry and is 
related to the free energy of isolated static fundamental charge. 

1.2 Monopole condensation 

A magnetically charged operator is constructed in such a way that it adds a monopole field to 
a given configuration. Its expectation value is given by 

(M> = f d-5) 

where Z is the partition function for the action in presence of a monopole. To better cope with 
fluctuations, one instead calculates the quantity 

P = -^Mv) = (S)s-(S) s (1.6) 

in which (S) § is the average of the action with a monopole insertion weighted with the modified 
action itself. Therefore, two simulations have to be run for each value of p\ The parameter (ju) 
should be different from zero in the confined phase where magnetic charges condense and drop to 
zero at deconfinement where magnetic symmetry is restored. This drop corresponds to a negative 
peak of p. In the vicinity of the critical temperature, j8 ~ j8 c , p is expected to scale as 

pL-'/^/fi'/v^.^)) (L7) 
The critical exponent v should be equal to 1/3 for a first order transition. 

2. Simulation details 

We have simulated two flavours of adjoint staggered fermions using the RHMC algorithm 
Trajectories had a length of NmdSi = 1, and an integration step of 8t = 0.02. Acceptance 

was above 90%. Inversions of the fermionic matrix were performed with the Conjugate Gradient 

algorithm. 

We have run simulations with two different lattice sizes, i.e. L 3 , x L t = 12 3 x 4, 16 3 x 4. We 
have calculated the chiral condensate, the Polyakov loop and the average plaquette for 12 values 
of j6 in the range (3.0,7.0). We computed the average plaquette and the p parameter. For time 
constraints, we were limited to runs at a bare quark mass of am = 0.04. In order to simulate an 
action with a monopole contribution, C* boundary conditions had to be implemented [0]. 

Our code has been run on the ApeMille machine in Pisa and the ApeNEXT facility in Rome. 

3. Results 

For each value of /3, we have calculated the Polyakov loop and the p parameter on both 
volumes, 12 x 4, 16 3 x 4. For the Polyakov loop, we reproduced the results of [|l|] at the same value 
of the bare quark mass, am = 0.04, see Figure (Q). Figure (|l|) shows the results for the p parameter. 
The position of the negative peak shifts slightly with the volume as one would expect for a physical 
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Figure 1: The p parameter, with am = 0.04, L f = 4, for two different spatial volumes 



transition. This estimate of the pseudocritical coupling is consistent with the position of the jump 
in the Polyakov loop. Furthermore, as j3 — > 0, one expects p to become volume independent, as 
H should tend to a nonzero constant. Our results are compatible with this expectation. We are 
currently working on the analysis of the scaling properties of p. We could not study the behaviour 
of p at the chiral transition as the bare quark mass turned out to be too large to see any visible effect 
in the chiral condensate susceptibility. 



4. Conclusions 



We have studied monopole condensation in lattice QCD with Nf = 2 staggered fermions in 
the adjoint representation. We evaluated the expectation value of a magnetically charged operator, 
which is expected to be an order parameter for confinement. For this operator, we observed a 
qualitative agreement with the expected behaviour. In particular, the data were consistent with 
the theoretical predictions for the /3 — > limit: p becomes volume independent. We are currently 
working on simulations with lighter quark masses and larger spatial volumes, in order to study the 
behaviour of the p parameter at the chiral transition and its scaling properties. 
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Figure 2: The Polyakov loop, with am — 0.04, L t — 4, for two different spatial volumes. Lines are left to 
guide the eye. 
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